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cd . Abstract 

Based on the linearized Eilenberger equations, the upper critical field {Hc2) of 



o 

O ' mixed d- and s-wave superconductors has been microscopically studied with 



an emphasis on the competing effects of mass anisotropy and spin Zeeman 

> 

CN ' coupling. We find the mass anisotropy always enhance Hc2 while the Zeeman 

o 



interaction suppresses Hc2- As required by the thermodynamics, we find Hc2 



Q\ ' is saturated at zero temperature. We compare the theoretical calculations 

with recent experimental data of YBa2Cu307_5. 
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It has been a consensustJ that a high-Tc cuprate superconductor has a (i-wave pairing 
symmetry, and CUO2 plane is responsible for superconductivity. However, a pure rf-wave 
symmetry is appropriate only for a tetragonal lattice structure, an orthorhombic material 
such as YBa2Cu307_5 (YBCO) is believed to have a subdominant s-wave component in 
the order parameter. The orthorhombicity in YBCO is originated from a mass anisotropy 
(MA) along a- and 6-directions. In other words, ma is larger than m},. Such a discrepancy 
is mainly due to CuO chains in the 6-direction.0 

Based on the Ginzburg-Landau theory, Xu et alia have explored effects of the mass 
anisotropy to show that an s-wave component always coexists with a dominant d-wave 
component in the bulk order parameter. Belzig et a/.Q have obtained the phase diagram of a 
mixed d- and s-wave superconductor in the quasiclassical theory,H and shown that the mass 
anisotropy gives rise to a non-zero s-wave component. 

In this paper, we shall investigate effects of the MA on the upper critical field (-^02) of 
a mixed d- and s-wave superconductor, based on the quasiclassical theory, using a quantum 
mechanical method we have developed in Ref.Q and compare a theoretical result with recent 
experimental data.Q We also take into account the paramagnetic Zeeman interaction (ZI) 
because Hc2 of a high-Tc superconductor is large at low temperature. As in Ref.Ha, we 
neglect any other effects associated with the chains. For example, we assume that lattice 
constants along a- and 6-directions have the same value. 

Since the calculation of ifc2 is a quantum mechanical problem of a charged particle 
in a constant magnetic fieldj^ first of all we need to check if solvable is a problem of a 
charged particle with two different effective masses along x- and y-directions, irix and iriy, 
respectively, in a constant magnetic field H = V x A = Hti. In the symmetric gauge, the 
Hamiltonian Ti of the particle is given by 

n= ^(Px + -Hyf + ^{py- -Hxf 
2mx^^ 2 ^' 2m/-^^ 2 ' 

1 eH e^H'^ 

(pI + ApJ) + T^iyp. - Xxpy) + (y^ + Ax^), (1) 



where A = rrix/my. Introducing operators a^ and ay such that 



and 
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'm^ujoVXx + ip^), (2) 



y/2m^Uo 
where ujq = \e\H/2mx, Ti. becomes 



S = /n^ . i^^'^oy + i'^Py), (3) 



H = u}o[a^ax + a^ay + V^ + i{a^ay - a^Ox-)]- (4) 



Let us introduce a new operator b = (a^ + iay) / y 2\^ to simplify H; then, we obtain 

n = iu,(bn + ^y (5) 

with [b,b'^] = 1, [b,7{\ = Ucb, and [&^,7-^] = —Ucb'^, where uj^ = \e\H / yffrQfn^. It shows 
that the problem we are considering reduces to that of a simple harmonic oscillator, and 
consequently Hc2 for an anisotropic mixed d-and s-wave superconductor can be exactly 
calculated. 

In the calculation of ifc2, we choose [110] and [—110] as x and y axis, respectively. An 
electron spectrum e(k) then changes from kl/2ma + kl/2mh to k"^ /2m + ck^ky/m. Here we 
have defined rUa = m/(l — c) and rrife = m/(l + c), where a small quantity c is introduced to 
represent the degree of mass anisotropy. The Fermi surface (FS) now is elliptical, and the 
density of state N^cj)) on FS becomes A^o(0)/[l + csin(20)]. By considering the anisotropy 
effect as a perturbation, we may assume that the order parameter can be still expanded 
in terms of the set of unperturbed eigenstates {/Ar(R)}. In Ref.Q, we have shown that the 
(singlet) order parameter is written as Aq 10)^+^2 12)^ + ^4 14)^ near Hc2, where we denote 
\N)d{s) as eigenstates f^ of d{s) channel. However, we expect that some other states such 
as |0)s, \l)s{d), and \2)d may involve in the order parameter due to the perturbation. The 
linear combination of these states for the order parameter near Hc2 will be determined by 
the symmetry of the system. 

The linearized Eilenberger equation is written as:El 



L/(cu, k, R) + isgn(cu)/i*H ■ [o-/(cu, k, R) - f{u, k, R)(t''] = 27rA(R, k) , (6) 

with 

A(R, k) = T ^(\/(k, k')f{u;, k', R))j.s , (7) 

UJ 

and 

L = 2\uj\+isgn{uj)vp-Il. (8) 

Here /(u;, k, R) is a quasiclassical Green's function, cr = {(Tx, Cy, Cx) are Pauli matrices, u is 
the Matsubara frequency, ^q can be interpreted as an effective magnetic moment of a quasi- 
electron with a mass anisotropy {rria 7^ iTih), which will be considered as a phenomenological 
parameter associated with a coupling strength between an electron spin and a magnetic 
field. The symbol {■ ■ ■)sf = I -^N^cj)') ■ ■ ■ represents the angular average over the Fermi 
surface. Since wp = vpa^ + vpb^ and 11 = — iVR — 2eA = Ilaa + Ilfeb with a and b being 
unit vectors along the a and h directions, 

Vi? ■ n = {vfx + CVFy)^x + {Vpy + CVFx)^y , (9) 

with VFx,y = kpx^y/rn in the x-y coordinate system. The pairing interaction in this coordinate 
system can be written as 

K(0,0') = Vs + Vdsin(20) sin(20O, (10) 

where = tan~^(A;j^/A;^). 

For the singlet pairing, / = f^iay and A = AoicTy. Using the inverse of the operator 
Lop = L + isgn{uj)^QH, which admits the representation 

/■oo 

Lop = / dsexp{-sLop) , (11) 



we show 



/•oo 
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Substituting Eq. (|T^) into Eq. (|^), we obtain the linearized gap equation of an anisotropic 
mixed d- and s-wave superconductor as follows: 

Ao(R,0) = 27rT5](y(0,0') /'rf^e-^[2|c.[+isgn^v^.n] 

xcos(2/i*iJOAo(R,0'))5F. (13) 

It is easy to see that Ao(R, 4>) turns out to be As(R) + Ad(R) sin(20) because of the pairing 
interaction V(0, 0') in Eq. (^). 

Let us, first of all, consider equations to determine Tc of such a superconductor. Setting 
H = 0, we obtain 

A. = N„V. (l + ^) lu Hi^A, - N,V. (£) m Hi|£A.. (14) 

and 

In order to calculate Tc of an anisotropic mixed d- and s-wave superconductor, it is convenient 
to introduce the transition temperature Td[s) of the anisotropic (i(s)-wave superconductor 
such that NoVs (l + y) In ^|^ = 1, and NoVd (| + fc^) In ^^ = 1. Then, T, can be 
expressed in terms of Tii{s) as follows: 

where c = cX]n=o(^ + V^)^^ with cj/j = {2n£) + l)7rTc. Here we would like to point out that 



a phenomenological value T^ ^ |(T(i + T^) + ^J [Td — TsY + Sc^T^Tg can be achieved only if 
Tg c:=LTd — Tc- If Ts << Trf, Tc has to be numerically calculated. 

Following Ref.H, we expand Arf(5)(R) in terms of {/Ar(R)}, namely, Ad^s(R) = 
Eiv^!^ /7v(R), to obtain 

N 



A,(R)= KE^ff I diY.Y.^nAi.N)c^lnii,c)f^.^+^{:R) 

N m n=0 

+KE A'7 ^eEE'^",™(e,iv)c^S(e,c)/^-n+,n(R) , (17) 



N rn n=0 



and 



A,(R)= Va T^^N fd^T^T. ^n,m{t N)C^nt (^' c)/^_„+„(R) 

N m n=0 

r ^ 



(18) 
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where 
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<^„,™(e,iV) = 27rTiVo5:e-2H«e " Vv^^ cos(2/iS/70 
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\e\H{vFi) 



2]{n+m)/2 
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if (n + m) is even and $n,m(C, N) = if (n + m) is odd, and 

C^^UiO- /^</> sin(20)^- [ ..^(vfC 

^nA^^^)-J 2vr[l + csin(20)]'^P\ ^^ [V2, 



[2csin(20) + c^ 



with C^% = Ci^%, C^^l^ = Ci'^^ and C^^J^ = Ci%. (See Appendix A for the detailed 
derivation.) It is necessary to investigate the symmetry properties of C^-'^ to calculate Hc2 
as we have mentioned early. Since c is a small quantity, we may expand the integrand of 
Cnm up to the c^ order. After a careful investigation of C^}^, we find 

sin[(m — n)n] 



C(o) oc 

"'"^ {n — m, + 4){n — rn + 2){n — 'm){n — rn — 2){n — rn — 4) 



ci'l 



C(o) 



OC 



(n — -m + 6) (n — -m — 6) 



(19) 



(20) 



and 



"'"" (n-m + 8)(n-m-8)' 



(21) 



As one can easily see, C^}^ ^ only if \n — m\ = 0, 2, 4, Cj^}^ ^ only if |n — m| = 0, 2, 4, 6, 
and C^}^ 7^ only if \n — m\ = 0, 2, 4, 6, 8. We also find other symmetry properties of C^}^ 
~^mn if \n — m\ = 2, 6, and C^}^ = C^^„ if \n — m\ = 4, 8. These properties 



such as Clj'l^ 



are still valid even if we expand C^}^ up to the c'^ order. We, thus, know that \2N)ii(s) do 
not couple to |2A^ + l)d{s) but to \2N)fi{s), and vice versa. Consequently, we expect that 
the order parameter is represented by J2N[^2N\'^^)d sin(20) + A2^\2N)s] near Hc2', in other 



words, Hc2 is determined by such coefficients as Agjy and A2iy. 



Using the ortho normality of {Jn}, we obtain the equation for the coefficients A)^ and 



(^) 



A]y from which Hc2 will be calculated as follows: 



N 



N -^ n=0 

r ^ 



(0) 
m+n—N 



(1) 
m+n—Ni 



(22) 



and 



n=0 



N 



N •' 71=0 

f ^ 



(1) 
m+n~N 



(2) 
m+n~N 



(23) 



71=0 



(1) 

riym+n—Ni 



with m + n - iV > 0, and |m - A^l = 0, 2, 4 for Ci°i,+„_jv, \m-N\= 0, 2, 4, 6 for C, 

and \m — N\ = 0,2,4,6,8 for C!^m+n-N- The equations to determine Hc2 are given by 

m = 0, 2, 4, ■ ■ -, namely. 



Ai'^=VsJd^ 






-|-/iO '±'o,0'-^0,0 "T ^2 ^2,0'-'2,0 "T ^4 ^4,o'-"4,0 



A['^=Vsjdi 



A^o^^lAi + ^(^^ E *Mcg^ + Af ) E ^Ufi% 



j=0 



i=0 



(srf) 
i,i 



i=0 



+<^K2cl,f + Af E ^Mf + 4'^ E *2+M^2: 

i=0 

Ai'^=V,Jd^ 



2 4 

^0 "^0,4^^0,4 + ^2 Z^ '^j,2+jQ,2+J + ^4 Z^ "^j.jQ,* 



i=0 j=0 

2 4 



i=0 



and similar equations of Ag , Ag , ^4 , ■ ■ -. Hc2 is the largest value of solutions which satisfy 
the condition for a non-trivial solution to exist in these equations of Ayv'* • ^^ ^^ have 
mentioned, \2N) (i{s){N = 0, 1, 2, • ■ ■) are involved in the determination of Hc2'-, however, it is 
expectable that the first few states such as \Q)d{s)) |2)d(s) and \'^)d{s) are important because 
in the case of c = 0, |0)d, \2)s and 4)^ play the dominant role in determining Hc2- Inclusion 
of more states such as \Gi)d(s) and \^)d(s) gives rise to a difference much less than 1%. 

As we did in Ref.Q, we introduce dimensionless unit in the calculation of Hc2'- t = 
T/Tc and h = 2\e\H {vp /27TTcy . In addition to the anisotropy parameter c, we also define 
parameters 6 = T^/Td and 7^ = {27ifiQ/evp)Tc, which is the strength of spin-magnetic 
field coupling. For the sake of comparison with experiment measurement, we convert the 
normalized magnetic field into the dimensional one by using 

g(^) = m. (24) 

_dH_\ rp _dh\ ■ V^^// 

By solving the eigen-equations for ^4^''' , we plot Hc2 in Fig. 1 as a function of temperature 
for several typical cases. Hc2 for an anisotropic mixed d- and s-wave superconductor with ZI 
taken into account is plotted with the solid line. Hc2 for an anisotropic mixed d- and s-wave 
superconductor without ZI is represented by the dot-dashed line. The involved parameter 
values are taken to be: 6 = Tg/Td = 0.06, c = 0.16, and 7^ = 0.15. Also plotted are Hc2 for 
a two-dimensional s-wave supercondcutor (dotted line) and a pure (i-wave superconductor 
(dashed line) without the MA and the ZI. Note that, even though the increase of 6 can 
enhance Hc2P for the value of 5 chosen here, the upper critical field for an isotropic mixed 
d- and s-wave order parameter (with 6 = 0.06) without ZI is more or less same as that 
of the pure (i-wave superconductor (dashed line). As shown in Fig. 1, we find that the 
mass anisotropy enhances Hc2 while the spin Zeeman coupling suppresses Hc2- This means 
that the mass anisotropy supports superconductivity; in other words, it increases T^ (and 
consequently the gap), which can be easily seen from the equation for Tc. In addition, one 
can also see that Hc2 is saturated at zero temperature as a reflection of the thermodynamic 
requirement; namely, on the phase boundary in the T ~ H plane, 

8 



dH,2 SS 

at T = near Hc2-, where 5S is the entropy difference between the normal and supercon- 
ducting states, and 5M is the magnetization near Hc2- 

Recently, the upper critical field Hc2 (parallel to the c axis) of YBCO with T^ = 84.3 K 
has been measured down to about AK.u The corresponding data is shown by solid circles 
in Fig. 1. The slope dH/dT\x^ has been experimentally observed to be —1.9 T/K.a In the 
theoretical calculation, we are mainly concerned about the low-temperature data of Hc2 most 
because at the low temperature thermal fluctuation effect is negligible. As shown in Fig. 1, 
the experimental data can be fit very well, with the above given parameter values, by our 
calculation for an anisotropic mixed d- and s-wave superconductor including the ZI. Here 
we would like to point out that, since 5 enhances Hc2 more significantly than cp if we take a 
large value of 5, we have to choose a physically unacceptable high value of 72 to fit the low- 
temperature data. Actually, it is believed that T^jEp ~ 0.1 for a high-Tc superconductor ,113 
which yields 7^ = t^Tc/2Ef ~ 0.157 by assuming /ig = hb (Bohr magneton). On the other 
hand, most of available experiments seem to indicate that 5 should be very small. All this 
facts demonstrate that the chosen set of parameter values are physically reasonable. Finally, 
we would like to mention: (i) Because Hp for the sample is about 185 T, the experimental 
data in Ref.Q are within paramagnetic limitlill with the critical ratio Hp/T^ estimated to be 
2.2 T/K for a rf-wave superconductor. Therefore, in our consideration, we do not have to 
include the spin-orbit interaction because it reduces the pair-breaking effect of the Zeeman 
interaction, and consequently, it allows ifc2 to be larger than Hp!^ In the theoretical point of 
view, since the strength of the spin-orbit coupling is proportional to Z^, where Z is an atomic 
number, we may neglect its effect in YBCO as long as no heavy-atomic impurity is taken 
into account as in this paper. However, it may play an important role in such heavy-fermion 
superconductors as UBeis and UPts.EJ (ii) Magnetic and non- magnetic impurities are pair 
breakers so that it is clear the impurities reduce ifc2 as well as T^.tJ However, the impurity 
concentration in the sample prepared in Ref.l3 seems to be negligible, the corresponding effect 



is not considered here, (iii) Recently, O'brien et al.iB have interpreted the experiment based 
on a three dimensional s-wave model.tJ However, it is well-known that superconductivity 
in YBCO is of two dimensional nature.^ (vi) A small deviation of theoretical results and 
experimental data occurs in high temperature region because thermodjTiamic fluctuations^^ 
of vortices is strong when temperature is high. 

In summary, the upper critical field (i^c2) of a mixed d- and s-wave superconductor with 
a mass anisotropy has been microscopically calculated based on the quasiclassical theory. We 
found the mass anisotropy supports Hc2 against Zeeman suppression. Hc2 becomes saturated 
at zero temperature in consistence with a thermodynamic requirement. The theoretical 
results are compared well with recent experimental data of YBa2Cu307_5. 

One of us (W.K.) would like to thank V. Kogan and S. C. Lee for helpful discussions. 
This work is supported by Texas Center for Superconductivity at the University of Houston 
and by the Robert A. Welch Foundation. 

APPENDIX A: 

Since Ao(R, 0) = As(R) + Arf(R, 0), we have two equations; namely, one is for As(R) 
and the other for Af;(R) as follows: 

A. W= y.EJ fn ii^cLm)^ I '^'~'"" -(2^;«0«(n.. n-) A.(R) 



27r[l + csin(20)] 
d(j) sin(20) 
2^[l + csin(20)] 



^"^ S / g [i;:itLl / ''^^-""" -(2.;^0"(n.. n-) A.(R) (Ai) 



and 



A.(R)- y^^l g |i;:ttl)i / ''^^-*" ^os(2KHom.,n^)^m 



^"^ ? / 'ijrfSw]I '^''""" ^os(2KHmin., n-)A.,(R), (A2) 

where 

fi(H+,H_) = 27rTA^oexp[-isgn(cu) fe] ^{{e-'^ - ice'^)U+ + {e''^ + ice-''^)U_}] (A3) 
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with n± = (n^. ± my)/^/2. Note that 



Hsgn(^) te j e(e-^* - zce^*)n+, -zsgn(a7) fe j ^(e^* + ^ce-^^)n_ 



2|e|i7 






[l + 2csin(2(/)) + c^ 



(A4) 



then we obtain 



a,(r)=kE 

a; 

+KE 



27r[l + csin(20)] 
d(j) sin(20) 
2^[l + csin(20)] 



d^e-^l"^'^ cos(2/i*i70^(n+, n_)A,(R) 
rf^e-^l"^!^ cos(2/i*ffO^(n+, n_) Arf(R) 



and 



(A5) 



A,{R)=V,Y: I ^TT^^^^T^ /rfee-^l-l«cos(2^*i70f^(n+,n_)A.(R) 



+v,E 



27r[l + csin(20)] 
d(f) sin(20)2 
2^[l + csin(20)] 



rf^e-'l^l^ cos(2/i*/70^(n+, n_)Ad(R), (A6) 



where 



fi(n+,n_) =27rriVoexp 



\e\H 






;i + 2csin(20) + c2 



xe 



-iSgn{Lj)(^^y{e-''l'-ice''t-)n+ -iSgn{<^)(^^)^(e'<^+ice-"^)n_ 



(A7) 



Expanding A5(d)(R) in terms of {/Ar(R)} and noting 



UA 



v"-/ 



h = j2\e\H 



^ VnTt\ ( In+i^ 



we obtain Eqs. (O) and (IT 



V 



N 



V ^^-1 J 



(A8) 
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T[K] 
FIG. 1. The upper critical field Hc2 for an anisotropic mixed d- and s-wave superconductor 
(solid line) with the spin Zeeman interaction taken into account. Solid circles are experimental 
data of YBCOQ. Also plotted are Hc2 of a two-dimensional isotropic s-wave superconductor (dotted 
line), a pure isotropic d-wave superconductor (dashed line), and an anisotropic mixed d- and s-wave 
superconductor (dot-dashed line), all obtained without the Zeeman interaction included. 
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